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In this paper a first order analytical system of difference equations is considered. 
For an asymptotically stable fixed point x° of the system a gradual approximation 
of the domain of attraction (DA) is presented in the case when the matrix of the 
linearized system in x° is a contraction. This technique is based on the gradual 
extension of the "embryo" of an analytic function of several variables. The analytic 
function is a Lyapunov function whose natural domain of analyticity is the DA and 
which satisfies an iterative functional equation. The equation permits to establish 
an "embryo" of the Lyapunov function and a first approximation of the DA. The 
"embryo" is used for the determination of a new "embryo" and a new part of the 
DA. In this way, computing new "embryos" and new domains, the DA is gradually 
approximated. Numerical examples are given for polynomial systems. 
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1 Introduction 



We consider the system of difference equations 



where g : — * f2 is an analytic function denned on a domain f2 included in 
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A point x° E f2 is a fixed point for the system (1) if x° satisfies 

x° = g(x°) (2) 

The fixed point x° of (1) is "stable" provided that given any ball B(x°,e) = 
{x E Q/\\x — x°\\ < e}, there is a ball B(x°, 5) = {x E Q/\\x — x°\\ < 5} such 
that if x E B(x°, 5) then g k {x) E B(x°, e), for k = 0, 1, 2, ... [1]. 

If in addition there is a ball B(x°,r) such that — > a; as fc — > oo for all 
x G B(x°,r) then the fixed point x° is "asymptotically stable". [1]. 

The domain of attraction Di(i ) of the asymptotically stable fixed point x° 
is the set of initial states x E Q from which the system converges to the fixed 
point itself i.e. 

DA(x°) = {xE n\g k (x) k ^ x } (3) 

It is known that a; is a fixed point for system (1) if and only if E M n is a 
fixed point for the system 

y k+1 = f(y k ) k = 0,1,2... (4) 
where / : Q — x° — > Q — x° is the analytic function defined by 

f(y) = g(y + x°)-x for y E Q - x° (5) 

The fixed point x° of (1) is asymptotically stable if and only if the fixed point 
E R™ of the system (4) is asymptotically stable. 

The domain of attraction of x°, DA(x°) is related to the domain of attraction 
of 0, DA(0) by the equation 

DA(x°) = DA(0) + x° (6) 

For the above reason in the followings instead of the system (1) we will consider 
the system (4). 

Theoretical research shows that the DA(0) and its boundary are complicated 
sets [2] , [3] , [4] , [5] , [6] . In most cases, they do not admit an explicit elementary 
representation. For this reason, different procedures are used for the approx- 
imation of the DA(0) with domain having a simpler shape. For example, in 
the case of the theorem 4.20 pg 170 [1] the domain which approximates the 
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DA(0) is defined by a Lyapunov function V built with the matrix dof of the 
linearized system in 0. In this paper, we present a technique for the construc- 
tion of a Lyapnov function V in the case when the matrix dof is a contraction, 
i.e. ||<9 /|| < 1. The Lyapunov function V is built using the whole nonlinear 
system, not only the matrix dof. V is defined on the whole DA(0), and more, 
the DA(0) is the natural domain of analyticity of V. The formula which de- 
fines the Lyapunov function V is used for determining an "embryo" of V and 
a first approximation of DA(0). The "embryo" is used for the determination 
of a new "embryo" and a new part of DA(0). In this way, computing new 
"embryos" and new domains the DA(0) is gradually approximated. 



2 Theoretical results 

Let be / : Q — > f2 an analytic function defined on a domain f2 C R n containing 
the origin € R n . 

Theorem 1 J/ t/ie funtion f satisfies the following conditions: 

/(0) = (7) 



H/H < 1 (8) 
then an asymptotically stable fixed point. DA(0) is an open subset of Q 
and coincides with the natural domain of analyticity of the unique solution V 
of the iterative first order functional equation 

V{f{x))-V{x) = -\\xf 
V(0) = 

The function V is positive on DA(0) and V(x) +oo, for any x° £ 

FrDA(0) (FrDA(0) denotes the boundary of DA(0)). 



PROOF. Let be a such that ||$o/|| < ck < 1. By the continuity of x h- > \\d x f\\ 
there exists a 5 > such that ||9 X /|| < a for ||x|| < 5. The mean value theorem 
gives 



\\f(x')-f(x")\\<a\\x'-x"\\ (10) 
for any x' and x" in the ball -5(0, S). Therefore 

\\f k (x)\\ <a k \\x\\ (11) 
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for any x in the ball -8(0, 5) and k — 0, 1, 2, .... For a ball -8(0, e) we take 5' = 
mm(e,S) and the ball B(0,5'). We have f k (x) G B(0,e) for any x G B(0,5') 
and k = 0,1, 2..., which means that is stable. 

From (11) we obtain f k (x) for any x G -8(0, 5). We can conclude now 
that is asymptotically stable. 

In order to show that DA{0) is an open subset of fl we consider x' from DA{0) 
and k > such that ||/ fco (£')ll < f • Because f k ° is a continuous function, there 
exists a ball B(x',S") such that \\f k °(x) - f k °(x')\\ < §, for any x G B(x',6"). 
Therefore, \\f k °(x)\\ < f < 5, for any x G 5"). It follows that x G DA(0) 

and therefore, B{x',5") C -DA(O). This proves that DA{0) is an open subset 

of a 

Now we consider x G ZL4(0) and the sequence {/ fc (x)}fc e p}. There exists k x such 
that / fe (:r) G B(0,5), for any fc > k x . Therefore, \\f kx+k (x) \\ < a k \\f kx (x)\\, 

oo 

for k = 0,1,2.... It follows that the series J2 \\f k { x )\\ 2 is convergent for any 
xeDA(0). 

Let be V = V(x) the function defined by 

oo 

V(x) = Y,\\f(*)\\ 2 iorxeDA(0) (12) 

k=0 

The above function defined on DA(0) is analytical, positive and satisfies (9). 
In order to show that the function V defined by (12) is the unique function 
which satisfies (9) we consider V = V'(x) satisfying (9) and we denote by V" 
the difference V" = V - V. It is easy to see that V"(f(x)) - V"(x) = 0, for 
any x G DA(0). Therefore, we have V"(x) = V"(f k (x)) for any x G DA(0) 
and any k = 0,1,2.... It follows that V"(x) = lim V"{f k {x)) = for any 

x G DA(0). In other words, V(x) = V'(x), for any x G DA(0), so V defined 
by (12) is the unique function which satisfies (9). 

In order to show that V(x) x - z ^> oo for any x° G FrDA(0) we consider x° G 
FrDA{0) and r > such that ||/ fe (a; )|| > r, for any k = 0,1,2.... For an 
arbitrary positive number JV > we consider the first natural number k\ 
which satisfies k\ > ^ + 1. Let be r\ > such that ||/ fc (x)|| > ^ for any 
k = 1,2, ..,ki and x G -E^rr ,?"!). For any x G B(x Q ,r\) fl DA{0) we have 

E ||/ fc (a:)|| 2 > N . Therefore, V(x) x -=^° oo. □ 

fc=0 

Remark 2 Newton's method for solving systems ofn nonlinear equations with 
n unknowns always leads to a system of difference equations which satisfies the 
conditions of Theorem 1. 
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We consider the expansion of V in 0: 



V(xi, x 2 , x n ) — ^2 E Bjij2-j n aJ i x 2 ''■■■ x3 n (13) 

m=2 \j\=m 



The domain of convergence of the series (13)is the set of those x G f2 which 
have the property that the series (13) is absolutely convergent in a neigh- 
borhood of x [7]. We denote by D° this domain. Actually, the domain of 
convergence D° coincides with the interior of the set of all points x° in which 
the series (13) is absolutely convergent [7]. 

Theorem 3 (Cauchy-Hadamard see [7]) A point x belongs to D° if and only 
if 

if JE \B~~~xfxfZx¥\<l (14) 
V b1=™ 

By the previous theorem we find a part D° of DA(0). 

In practice, we can compute the coefficients -Bj U2 ...jn U P to a finite degree p. 
This degree p has to be big enough for assuring that the domain D® given by 

D° p = {xe SI/ \B jlj2 ... jn x?x?...xt\ < 1} 

V \j\=p 

approximates the region of convergence D° of the series of V and the "embryo" 

v p °(x u x 2 ,..., Xn ) = E E i[rj,.,y:-' J 2---< (is) 

rn=2 \ j\=m 

approximates V with accuracy. 

The first estimate of the DA will be D®. 

In order to extend the first estimate D®, we expand V in a point x° close to 
FrDp in which | V^°(rc°) | is still small. That is because, according to Theorem 
1, the points x close to FrD® for which |V p °(x)| is extremely high are close to 
FrDA(0). 

To find the expansion of V in x° close to FrD®, we will compute the expansion 
in x° of the "embryo" VI of V. We obtain: 
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V p 1 (x l ,x 2 ,...,x n )=j: £ . ./ (x°)(x 1 -x^...(x n -x° n )^/ml 

=0 \j\=m aX ^ aX 1 ■■■ ax n 



m- 
V 



m=0 \ j\=m 



We consider the set 



Dl = {xe SI/ J Yl \ B ,m...jn(xi - x i)H x 2 - x ^...{x n - x°y»\ < 1} 

V \3\=P 



which provides a new part D p of DA(0). 

So, U Dp gives a larger estimate of DA(0). We can continue this procedure 
for a few steps, till the values \V p k \ become extremely large and we obtain the 
estimate D° p U D\ U ... U D k v of DA(0). 



3 Numerical results 



3. 1 Systems with known domains of attraction 



In this subsection, we will present some examples of systems of one, two or 
three difference equations, for which we can compute easily the DA. We will 
apply our technique to these examples, and we will show how the real do- 
mains of attraction are gradually approximated. These examples are meant to 
validate our procedure. 

The computations were made using a program written in Mathematica 4, 
Wolfram Research on an Intel Pentium III PC (2Ghz, 512MB of RAM). The 
data for the estimations (the degree up to which the approximation is made, 
the necessary timing for the estimations) are displayed in Table 1. In our 
figures, the thick black line represents the true boundary of the domain of 
attraction, the dark grey set denotes the first estimate Dp of DA and the 
further estimates D k of DA with k > 1 are colored in light grey. 
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Example 1 



We consider the following one dimensional difference equation: 



x n+l 




(17) 



This difference equation has two fixed points: x = and x = —0.271845. 
The fixed point x = is asymptotically stable, while the other fixed point 
is unstable. It can be proved (by the staircase method) that DA(0) is the 
interval (-0.271845,0.653564). 

We applied the procedure described above to obtain an estimation of the 



• After the first step, we obtained D° = (-0.27184,0.27184). 

• After the second step applied in x 1 = 0.2718, we obtained D 1 = (0.01345, 0.53015). 

• After the third step applied in x 2 = 0.5, we obtained D 2 = (0.38378, 0.61622). 

• After the forth step applied in x 3 = 0.61, we obtained D 3 = (0.59785, 0.622175). 

Therefore, the estimate of DA(0) obtained after four steps is (-0.27184, 0.622175). 



Example 2 

In [1], the following system of difference equations is considered: 



It is easy to see that (x, y) = (0, 0) is an asymptotically stable fixed point for 
the system (18), and the matrix of the linearized system in (0, 0) satisfies the 
conditions from Theorem 1. There are two more fixed points of system (18) 
(which are unstable) which are represented by the gray points in Fig 1.1-1.2. 
The DA(0) is the set: 



Using the technique described above, after 2 steps of gradual extension of the 
"embryo", we find the estimate of DA(0) presented in Figure 1.1. If we apply 
the second step for several points close to FrD°, we obtain the estimate of 
DA(0) presented in Figure 1.2. 



DA(0). 



< 



Xn+l — — 2/n(l — % 2 n — Dn) 

^ Dn+i = —x n {l — x 2 n — y 2 ) 



(18) 



DA{0) = {(x, y) e M 2 /x 2 + y 2 < 3} 
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Fig. 1.1. The estimate of DA(0) after 2 Fig. 1.2. The estimate of the DA(0) after 
steps for system (18) 2 steps, in several points close to FrD° 

Example 3 

The following decoupled system of difference equations is considered: 



*n+l-44 (ig) 
Vn+l = 9Vn 

It is clear that (0, 0) is an asymptotically stable fixed point for the system (19) 
and 

™<o) = (-i.i)x(-i,i) 

Applying the above procedure, after only one step, we obtain the whole DA(0). 





Fig. 2. The estimate of DA(0) after 1 steps for system (19) 
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3.2 Systems for which we don't know the domain of attraction 

In this subsection, some systems of difference equations are presented for which 
we don't know the DA(0). For these examples, we will apply the procedure 
presented above. 

Example 4 

We consider the difference equation: 

x n+1 = -x 2 n - 2x 3 n - 44 - &4 (20) 

This equation has only one fixed point, namely x = 0. We obtained 

• D° — (-0.442585,0.442585), after the first step. 

• D 1 = (-0.673088, -0.212082), after the second step applied in x 1 = -0.44258. 
The estimate of the DA(0) obtained after two steps is (-0.673088,0.442585). 

Example 5 

The following system is considered: 

Kn+l ^"n-Vn (21) 

This system of difference equations has besides the asymptotically stable fixed 
point (0,0), an unstable fixed point represented in the following figure by 
a black point. After two steps, we can obtain the following estimate of the 
DA(0): 
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Fig. 3. The estimate of DA(0) after 2 steps for system (21) 

Example 6 

The following system of three difference equations is considered: 



Vn+1 



1 1 19 19 1 9 1 9 1 

2 J - , nyn ~ ^n^n ~T~ 3 J 'n"ra ' 12 n n 3 w"n 12 i n 12 "nan^n 

1 1 12 1212 121 

1 1 12 12 12121 



[22) 



The estimate of the DA of the asymptotically stable fixed point (0, 0, 0) ob- 
tained after one step is presented in Fig. 4. 




Fig. 4. The estimate of DA(0) after 1 step for system (22) 
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Table 1. Numerical data 



example 


order of approximation p 


timing for I s * step 


timing for 2 nd step 


1 


4096 


9.3 h 


32.5 h 


2 


164 


10 min 


24.1 h 


3 


500 


12.4 min 




4 


625 


1.2 h 


12.4 h 


5 


256 


16.8 min 


45.2 h 


6 


54 


4.6 h 
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